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Abstract 

We extract atom-atom potentials from the effective spherical potentials that 
successfully model Hugoniot experiments on molecular fluids, e.g., O2 and iV~2. 
In the case of O2 the resulting potentials compare very well with the atom- 
atom potentials used in studies of solid-state properties, while for N2 they 
are considerably softer at short distances. Ground state (T = OK) and room 
temperature calculations performed with the new N — N potential resolve the 

previous discrepancy between experimental and theoretical results. 
PACS numbers: 31.70.K, 34.20. G, 61.25.E 
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In recent years there has been an increasing interest in the high-pressure properties of 
liquids and solids HJ. As experimental research focused on the 'extreme' states of matter, 
it also forced us to refine our knowledge of the relevant intermolecular interactions. In 
many instances, including the study of thermodynamic properties, such interactions can be 
described by effective two-body potentials. Because even simple molecular systems such as 
N2 or O2 can reveal bewildering complexity at high pressures, developing reliable potentials 
has proven to be a laborious task. A further difficulty is that the intermolecular interactions 
are in most cases anisotropic. Taking into account the anisotropy is essential for an accurate 
description of solid phases. In general, site-interaction models (commonly referred to as 
atom-atom potentials), describing the intermolecular interaction as a sum of the interactions 
between atoms residing on different molecules, have proven to be a successful approach. 

Quantum mechanical calculations have been used to provide information on the inter- 
molecular atom-atom potentials, which are further adjusted to reproduce measured quan- 
tities such as virial coefficients, sublimation energy, melting point at room temperature, 
low-pressure crystal structures, etc. These experimental data are mainly suitable to 

constrain the potential in the attractive region. There are also abundant high-temperature 
and high-pressure data which could provide information on the short-range anisotropic re- 
pulsion. However, they have been mostly interpreted assuming a spherical short-range re- 
pulsion, with addition of a long-range van der Waals attractive tail f§-§ . The main reason 
is that the statistical mechanical theory of particles interacting with spherical potentials 
is well developed, and reliable and efficient computational tools are available. In contrast, 
such theories for anisotropic interactions are cumbersome, requiring approximations that 
are hard to justify and generally limited in their applicability 0. 

The fact that the spherical short-range repulsions employed by earlier workers can suc- 
cessfully explain high-density, high-temperature experimental data of molecular fluids sug- 
gests that the sphericallization of the anisotropic repulsion is reasonable in the high-density, 
high-temperature environment. In this work we set out to investigate to what extent one 
can extract the anisotropic interactions, in particular at short range, from such spherical 
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interactions. Considering the extensive database of spherical intermolecular potentials for 
both like and unlike pairs, the atom-atom potentials so extracted could provide a unified 
description of pure and mixed molecular solids under a variety of pressure and temperature 
conditions 

Physically, the spherical potentials extracted from experimental data for anisotropic 
systems are effective interactions, resulting from 'averaging' the anisotropy at high temper- 
atures. A great deal of effort has gone into justifying and finding procedures to construct 
such effective spherical interactions. In this regard, the so-called 'median' potential of Shaw, 
Johnson, and Holian || represents one outstanding contribution. Subsequently, Lebowitz 
and Percus || clarified the main ideas behind this approximation, which we now summarize 
for the slightly more general case of a mixture. 

If y '(r, Qi, Q 2 ) is an anisotropic potential between molecules of type i and j, and 0o ( r ) 
is the corresponding effective spherical potential, a perturbation potential, 



V,fi 1 ,n 2 )=^(r) + A 7 [^,^], (1) 



is constructed to connect 0o( r ) at 7 = to 4>^(r, Oi, fi 2 ) at 7 = 1 with < 
7 < 1, A o [0^',0^] = 0, A 1 [0^',0o i ] = 4> ij ~ 4>o- The first-order Helmholtz 
free energy correction along this path is AF = y(dF / dy)\o and is proportional to 
J2ij XiXj j <7^-(ri, r 2 )(dA^ /d r y)\ dr 1 dr 2 dQ 1 dQ 2 , where g° 's are the pair correlation functions 
of the sphericalized system and Xi S £1X6 the concentrations (J2i x % = !)• ft is eas y t° see that 
AF can be annulled, independent of all g^s, by requiring 

J(dA^/d 1 )\ dQ 1 dQ 2 = (2) 

for all pairs. 

The disadvantage of Eq. (2) as a definition for spherical potentials is that (dA % j / dy)\o 
is, to a large extent, arbitrary. The median is defined by (dA^/dy)\o = sgn((fi l:j — (ffi ) 
and minimizes the integral of absolute potential deviations, J \(f>^ — <f)Q\dQidQ 2 |n|. We 



also note that, as a consequence of the Gibbs-Bogoliubov inequalities, for the common case 



A 7 [0* J , 0o ] — iW 1 — 0o ] the first order perturbation theory provides a rigorous upper bound 
on the free energy of the system of interest |7j. While for the general case that we consider 
here only approximate inequalities hold []10|1 , they do lead to useful heuristic variational 



principles. 

One of the main advantages of the median potential over other effective spherical poten- 



tials |Tl| is that it is independent of density and temperature (and also concentrations in 
the case of mixtures [0)1111 )• Despite its essentially heuristic basis PJIif, this appears to be 
an optimal choice for the thermodynamics of simple fluids such as N 2 and CO2 at high pres- 
sures and temperatures ||, where the repulsive interactions make dominant contributions. 
Hence, we treat the spherical potentials derived from experimental data (such as Hugoniot 
data for iV 2 and 2 ) as median potentials and develop a procedure to invert them to extract 
the corresponding anisotropic potentials. 

To illustrate the basic idea behind such an inversion process, consider first a simpler 
system of dumbbells consisting of two hard spheres (with diameter R) with their centers 
fixed at length L apart. In this case the equivalent problem is to determine R by inverting 
an effective hard-sphere potential (with diameter R m ). The median definition of R m is 
/ sgn(Ri2 — R m )dQ 1 dQ 2 = £L5|, where i?i 2 is the distance of closest approach between two 
dumbbells oriented at solid angles Qi and f2 2 . It is fairly simple to show that this equation 
has a unique solution R. 

Potentials that are used to model the thermodynamics of high-pressure, high-temperature 
fluids are more complicated than pure hard-cores. The most successful one is the so-called 
exponential-6 ('exp-6') potential, 

"W^j^Wi-^)]-"^)}, (3) 

that combines the short-range exponential repulsion found in early quantum-theory calcu- 



lations |16| with a truncated dispersive interaction [|17j that accounts for the attraction at 
longer distances. This functional form has been found to provide a good representation 
of the experimental Hugoniots of iV 2 , O2, and other small molecules [£§-|6|], while satisfying 



to a large extent the 'law of corresponding states' |4[]. In this work we use the potential 
parameters for nitrogen, e N2 jks = 101. 9K, r$ 2 = 4.09A, a N2 = 13.2, obtained by Ross 
and the oxygen parameters, e° 2 jks = 125K, r° 2 = 3.86A, a° 2 = 13.2, that have been also 
successfully used to model the molecular dissociation under shock compression As 
experimental Hugoniots are excellent probes of the short-range repulsive interactions, such 
potential parameters are perhaps the best representation of the repulsive region under the 
spherical approximation. We also note that the well position r$ and depth e of the above 
potentials are in very good agreement with the ones obtained in molecular beam scattering 



experiments [IS 



The present inversion process derives an anisotropic potential 0(r, Q 2 ) that maps into 
v(r) (with the exp-6 potential parameters described above) through the median construction, 

J sgn[(j){r, Q. x , fl 2 ) - v(r)]dfiit2fi 2 = 0, (4) 

at all intermolecular distances r, where <f>(r, Oi, fi 2 ) is the sum of interactions (pfcj) between 
atoms i and j residing on different molecules, 

0(r,Qi,Q 2 ) = XM^-), (5) 

ij 

We also use an exp-6 form to represent <p(r) in Eq. (5) and recast Eq. (4) as a x 2 minimiza- 
tion. 

To elaborate, for a given set of exp-6 parameters (e,r , a) of <£>(r) [thereby a given 
0(r, Qi, Q 2 )} we calculate the spherical potential 0°(r; e, a, r Q ) that satisfies the median con- 
dition, Eq. (4), at a large number of intermolecular separations (~ 100), between r min and 
a large cutoff separation (r max ) beyond which the potential can be neglected. In practice 
we choose r min = 2 A, [v(r min )/k B w 10 5 K] and r max = 7 A, which is sufficient to cover the 
relevant potential region up to very high densities and temperatures. We then minimize 
the function x 2 (e, r ,a) = Z)j[ t '( r i) — 0°( r i; e > r o ; a )] 2 /[ v ( r i) 2 + e *] to determine the exp-6 
parameters of <p(r), where is the well depth of v{r) JOJ. 

We carried out such minimizations for 2 and A^ 2 , with fixed gas-phase bond lengths of 
lo 2 — 1.20741 A and = 1.097685A. The resulting potential parameters are: 



• 0-0 potential: e/k B = 45.64K, a = 12.52, r = 3.480A; 

• N-N potential: e/k B = 34.42K, a = 12.59, r = 3.773A. 

In Fig. 1 we compare the — potential with a potential that is used to study the 
solid-state properties of molecular oxygen 0. The agreement is satisfactory, in particular, 
at short distances. In Fig. 2 we make a similar comparison between the N — N potential and 
the so-called Etters potential which is widely used for the study of molecular nitrogen. 
Fig. 2 also shows another potential ||2U|| , which was recently used for the study of small N 2 
clusters |ZT| . It is clear that the potential that we extracted by inverting the median is much 
softer at small separations than either of these two, the main difference being attributable 
to the exp-6 parameter a that mediates the stiffness of the repulsion. 

Solid nitrogen has been the subject of numerous experimental and theoretical studies in 
the last decades leading to a good understanding of its phase diagram at lower pressures 



P]22|--[25f . However, questions remain regarding its high- pressure behavior, i.e., above ~ 20 
to 30GPa. Some of these questions concern the crystal symmetry and thermodynamics of the 
high-pressure phases, and others the stability of molecular nitrogen with respect to atomic 
dissociation and its possible transition to a metallic state. But ultimately all these questions 



appear to be related to each other p3| , |26| . X-ray diffraction and Raman spectroscopy studies 
of solid molecular nitrogen have revealed an increasing complexity of the phase diagram 
with pressure, e.g., the e — N 2 phase, which is apparently a high-pressure distortion of the 
Pm2>n cubic lattice of 5 — iV 2 [BJ. Calculations using a site interaction model such as the 
Etters potential have successfully described the low pressures part of the phase diagram. 



However, they show sharp disagreement with the experimental data at high pressure |27 
which in addition to a high sensitivity of the high-pressure structures to the intermolecular 
potential pSfl , emphasizes the need for more reliable interactions. To test the N—N potential 



obtained in this work, we first calculate the OK-pressure of e — N 2 (rhombohedral R3c) by 
minimizing at fixed density the energy with respect to distortions of the Pm3n structure 
consistent with the R3c symmetry fl3"|,|29|j . Because experimental data are available only at 



room temperature, we correct the OK-pressure, p(p, 0), by adding the thermal contribution, 



?>^pk B T, at each molecular density p (7 is the Griineisen parameter) P?| . In addition, we 
performed isothermal molecular dynamics simulations at room temperature. 

Fig. 3 compares these results with the experimental data of p7| and the results obtained 



with the Etters potential |3U] . The overall agreement is remarkably good, given the fact that 
the N — N potential has not been adjusted at all and is a significant improvement over the 
Etters potential calculations. This demonstrates the crucial contribution of the short-range 
repulsion to close-packed high-pressure structures. Our OK energy minimizations do not 
show any significant relaxation of the R3c structure to a lower symmetry, higher-pressure 



R3c lattice [p9|] . This may be due to the fact that we are only considering structures 
with eight molecules per unit cell and also without inclusion of the quadrupole-quadrupole 
interaction. However, such a transformation is believed to occur at a still higher pressure 



and low temperature [31 



We further test the agreement of the experimental and theoretical results by looking 
at the isothermal bulk modulus B = V(dP / 8V)t- In order to extract B we fit both the 
experimental and theoretical results with a modified Birch equation, that has been previously 
used to model the bulk modulus of nitrogen at pressures up to about 2GPa |52J : 



P = Po + 3 -^(y 7 -y 5 )[i-ay 2 -i)] (6) 

where y = (Vo/V)^. The comparison shown in Fig. 4 yields very satisfactory agreement. 

Local-density-functional energy calculations |25| yield a reliable description of molecular 
nitrogen, but one that is limited to OK. In order to further asses the accuracy of our N — N 
potential we calculate the OK-energy as a function of density for a — N2. Fig. 5 shows again 
good agreement between these data and the corresponding local-density-functional results 
p5| . The potential obtained in the present work gives a slightly steeper isotherm at high 
density, producing good agreement with experimental pressures in Fig. 3. 

We conclude that the N — N and O — O potentials extracted from the median provide 
a good representation of the short-range repulsive region. Given the fact that classical 



molecular dynamics and Monte Carlo simulations are still the most accurate and efficient 
way of studying anisotropic molecular systems at finite temperatures, the need for reliable 
potentials for such simulations cannot be overstated. The present work has shown that the 
task of constructing such anisotropic interactions is greatly simplified by inverting a spherical 
potential based on experimental data. The resulting potential is 'transferable', in that it can 
self-consistently describe material properties of different solid phases as well as liquid phases. 
If desired, one can refine the attractive part of the potential with inclusion of additional 
multipolar interactions |§. Such additional refinements that can effectively capture the 
effect of the long range two and three-body dispersive interactions may be necessary to 
describe properties at low pressures and temperatures. However, at very high pressures the 
structure and thermodynamics of the fluid is dominated by excluded volume effects well 
represented by two body potentials ||,[33[ that, as we show here, are also 'transferable'. 

We thank Andrew K. McMahan for kindly providing the results published in [p5| |. This 
work was performed under the auspices of the U. S. Department of Energy by University of 
California Lawrence Livermore National Laboratory under Contract No. W-7405-Eng-48. 
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FIGURES 

FIG. 1. Repulsive region comparison between the present O — O potential (circles) and that of 
ref. || (solid line). 

FIG. 2. Repulsive region comparison between the present N — N potential (circles) and that 
of ref. M (solid line) and ref. [pCf] (dashed line). 



FIG. 3. Room temperature isotherm of iV^: molecular dynamics (solid line) and energy mini- 
mization (dot-dashed line) results obtained with the N — N potential of this work; Etters potential 
calculations from ref. |3(| (dotted line); experimental data of ref. |27|] (circles: e — N2, triangles: 
8-N 2 ). 

FIG. 4. Bulk modulus B of N2 extracted from molecular dynamics results obtained with the 



N — N potential of this work (solid line) and experimental data of ref. [27| (circles: e — N2, triangles: 
6-N 2 ). 

FIG. 5. OK energy calculations for a — N2 with the N — N potential of this work (solid line) 
and the local-density- functional results of ref. [25|. 
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